Mon. Not. R. Astron. Soc. 000, 000-000 (0000) 



Printed 25 January 2012 



(MN M£X style file v2.2) 



Triaxiality, principal axis orientation and non-thermal 
pressure in Abell 383 



Andrea Morandi 1 *, Marceau Limousin 2 ' 3 

1 Raymond and Beverly Sackler School of Physics and Astronomy, Tel Aviv University, Tel Aviv, 69978, Israel 
(N. Laboratoire d'Astrophysique de Marseille, Universite de Provence, CNRS, 38 rue Frederic Joliot-Curie, F-13388 Marseille Cedex 13, France 
3 Dark Cosmology Centre, Niels Bohr Institute, University of Copenhagen, Juliane Maries Vej 30, DK-2100 Copenhagen, Denmark 

25 January 2012 

ABSTRACT 

While clusters of galaxies are regarded as one of the most important cosmological 
probes, the conventional spherical modeling of the intracluster medium (ICM) and 
the dark matter (DM), and the assumption of strict hydrostatic equilibrium (i.e., the 
equilibrium gas pressure is provided entirely by thermal pressure) are very approximate 
at best. Extending our previous works, we developed further a method to reconstruct 
(— I \ for the first time the full three-dimensional structure (triaxial shape and principal 

O ■ axis orientation) of both dark matter and intracluster (IC) gas, and the level of non- 

thermal pressure of the IC gas. We outline an application of our method to the galaxy 
cluster Abell 383, taken as part of the CLASH multi-cycle treasury program, presenting 
results of a joint analysis of X-ray and strong lensing measurements. We find that 
the intermediate-major and minor-major axis ratios of the DM are 0.71 ± 0.10 and 
0.55 ±0.06, respectively, and the major axis of the DM halo is inclined with respect to 
the line of sight of21.1±10.1 deg. The level of non-thermal pressure has been evaluated 
to be about 10% of the total energy budget. We discuss the implications of our method 
for the viability of the CDM scenario, focusing on the concentration parameter C and 
the inner slope of the DM 7, since the cuspiness of dark-matter density profiles in the 
central regions is one of the critical tests of the cold dark matter (CDM) paradigm 
for structure formation: we measure 7 = 1.02 ± 0.06 on scales down to 25 Kpc, and 
C = 4.76 ± 0.51, values which are close to the predictions of the standard model, 
and providing further evidences that support the CDM scenario. Our method allows 
us to recover the three-dimensional physical properties of clusters in a bias-free way, 
overcoming the limitations of the standard spherical modelling and enhancing the use 
of clusters as more precise cosmological probes. 

Key words: cosmology: observations - galaxies: clusters: general - galaxies: clusters: 
individual (Abell 383) - gravitational lensing: strong - X-rays: galaxies: clusters 



1 INTRODUCTION 

Clusters of galaxies represent the largest virialized struc- 
tures in the present universe, formed at relatively late times 
and arisen in the hierarchical structure formation of the Uni- 
verse. As such, their mass function sensitively depends on 
the evolution of the large scale structure (LSS) and on the 
basic cosmological parameters, providing unique and inde- 
pendent tests of any viable cosmology and structure forma- 
tion scenario. 

Clusters are also an optimal place to test the predic- 
tions of cosmological simulations regarding the mass profile 
of dark halos. A fundamental prediction of N-body simula- 
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tions is that the logarithmic slope of th e DM 7 asymptotes t o 
a shallow powerlaw trend with 7 = 1 l|Navarro et al.l fl996), 
steepening at increasing radii. Furthermore, the degree of 
mass concentration should decline with increasing cluster 
mass because clusters that are more massive collapse later, 
when the cosmological background density is lower (e.g., 
iBullock et"aT1l200ll : lNeto et al.ll2007T ). Measurements of the 
concentration parameters and inner slope of the DM in clus- 
ters can cast light on the viability of the standard cosmo- 
logical framework consisting of a cosmological constant and 
cold dark matter (ACDM) with Gaussian initial conditions, 
by comparing the measured and the predicted physical pa- 
rameters. With this respect, recent works investigating mass 
distributions of individual galaxy clusters (e.g. Abell 1689) 
based on gravitational lensing analysis have shown poten- 



© 0000 RAS 



2 Morandi et al. 



tial inconsistencies between the predictions of the CDM 
scenario relating halo mass to concentration parameter, 
which seems to lie above the mass-concentration relation 



predicted by the standard ACDM mod el (Broadhur st et al 



20051; [Limousin et alj|2007al : lOguri et al"j|2009t IZitrin et al. 
2011al lbh. though other works (e.g. lOkabe et alJl20ld ) found 



a distribution of the concentration parameter in agreement 
with the theoretical predictions. Lensing bias is an issue here 
for clusters which are primarily selected by their lensing 
properties, where the major axis of a cluster may be aligned 
preferentially close to the line of sight boosting the pro- 
jected mass density obser ved under the assumption of stan- 
dard spherical model ling llGavazzil 120051 ; lOguri fe Blandfordl 
l200gh . For example, iMeneghetti et all l|2010l ) showed that, 
owing to the cluster triaxiality and to the orientation bias 
that affects the strong lensing cluster population, we should 
expect to measure substantially biased up concentration pa- 
rameters with respect to the theoretical expectations. In- 
deed, departures from spherical geometry play an important 
role in the determina tion of the desired physical parameters 
<|Morandi et al.ll2010l ), and therefore to assess the viability 
of the standard cosmological framework. 

In particular, the conventional spherical modeling of the 
intracluster medium and the dark matter is very approxi- 
mate at best. Numerical simulations predict that observed 
alaxy clusters are triaxial and not spherically symmetric 



.Shaw et al.ll2006l;lBett et alj2007l : lGottlober fc Yepesl2007l ; 
Wang fc White! 120091 1. as customarily assumed in mass de- 
terminations based on X-ray and lensing observations, while 
violations of the strict hydrostatic equilibrium (HE) for X- 
ray data (i.e. intracluster gas pressure is provided entirely by 
thermal motions) can sy stematically bias low the determi- 
nation of cluster masses (lLau et al.ll2q09l;IZhang et al.ll2010l ; 
iRichard et"al]|2010l : IMeneghetti et alj|20ich . ~ 

The structure of clusters is sensitive to the mass density 
in the universe, so the knowledge of their intrinsic shape has 
fundamental implications in discriminating between differ- 
ent cosmol ogical models and to constrain cosmic structure 
formation (jlnagaki et al.lll995l ; ICorless et afll2009h . There- 
fore, the use of clusters as a cosmological probe hinges on 
our ability to accurately determine both their masses and 
three-dimensional structure. 

By means of a joint X-ray and lensing analysis, 
Mor andi et al. I (|2007l , l2010l . I 201 la, b) overcame the limita- 
tion of the standard spherical modelling and strict HE as- 
sumption, in order to infer the desired three-dimensional 
shape and physical properties of galaxy clusters in a bias- 
free way. 

In the present work, we outline an application of our 
updated method to the galaxy cluster Abell 383. This is a 
cool-core galaxy cluster at z = 0.189 which appears to be 
very relaxed. While in the previous works we assumed that 
the triaxial ellipsoid (oblate or prolate) is oriented along the 
line of sight, here we developed further our method by recov- 
ering the full triaxiality of both DM and ICM, i.e. ellipsoidal 
shape and principal axis orientation. We discuss the impli- 
cations of our method for the viability of the CDM scenario, 
focusing on the concentration parameter and inner slope of 
the DM. 

Hereafter we assume the flat A.CDM model, with mat- 
ter density parameter Q. m = 0.3, cosmological constant 
density parameter Qa = 0.7, and Hubble constant Ho = 



70 km s 1 Mpc 1 . Unless otherwise stated, quoted errors are 
at the 68.3% confidence level. 



2 STRONG LENSING MODELING 
2.1 Multiple Images 

Abell 3 83 have been prev i ously studied using lensing tech- 



niques JSmith et al 
l2010l; iNewman et al.l 



2001 



2011 



Sand et al l |2008|; IRichard et al 
Richard et all 1201 ll ; IZitrin et al 



2011b! ') . We benefit from these earlier works in order to pro- 
pose multiply imaged system, in particula r the most re- 
cent works based on the CLASH HST data (jPostman et al.l 
20111). We agree with the identifications proposed by 



Zitrin et al.l (|2011b! ) and we did not find any additional mul- 



tiply imaged system, although we notice some blue features 
in the core of the cluster that may be multiply imaged. We 
use 9 multiply imaged systems, five of them have a spec- 
troscopic redshift measurement. For the remaining systems, 
the redshifts will be let free during the optimization. They 
are listed in Table [1] 

2.2 Mass Distribution 

The model of the cluster mass distribution comprises three 
mass components described using a d ual Pseudo Isothermal 
Elliptical Mass Distri bution (dPIE iLimousin et al. I 120051 ; 
lEhasdottir et alj|2007l ). parametrized by a fiducial velocity 
dispersion a, a core radius r corc and a scale radius rj^: (i) 
a cluster scale dark matter halo; (ii) the stellar mass in the 
BCG; (iii) the cluster galaxies repr esenting local perturba - 
tion. As in earlier works (see, e.g. ILimousin et al J 12007a! ) . 
empirical relations (without any scatter) are used to relate 
their dynamical dPIE parameters (central velocity disper- 
sion and scale radius) to their luminosity (the core radius 
being set to a vanishing value, 0.05 kpc), whereas all geo- 
metrical parameters (centre, ellipticity and position angle) 
are set to the values measured from the light distribution. 
Being close to multiple images, two cluster galaxies will be 
mode lled individually, namely PI and P2 (jNewman et al.l 
l201ll ). Their scale radius and velocity dispersion are opti- 
mized individually. We allow the velocity dispersion of clus- 
ter galaxies to vary between 100 and 250 kms -1 , whereas 
the scale radius was forced to be less than 70 kpc in or- 
der to ac count for tidal stripping of their dark matter halos 
(see, e.g. ILimousin et a"ill2007bl . 120091 : iNataraian et al.ll2009l : 
IWetzel fc Whitell2010l . and references therein.) 

Concerning the cluster scale dark matter halo, we set 
its scale radius to 1 000 Kpc since we do not have data to 
constrain this parameter. 

The optimization is p erformed in the im age plane, using 
the LenstooiQ software (| Julio et al.ll2007l ). 



2.3 Results 

Results of the strong lensing analysis are given in Table [2] 
The RMS in the image plane is equal to 0.46" . In good agree- 



Hereafter we use the notation r s to indicate the scale radius in 
spherical symmetry, and R B for the triaxial symmetry. 
2 http:/ /www. oamp.fr/cosmology/lenstool/ 
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Figure 1. Core of Abell 383. Size of the field equals 55x55 arcsec 2 corresponding to 174x174 kpc 2 . Multiply imaged systems used 
this work arc labelled. 



ID 




R.A. 




Decl. 


■^spec 


^•modcl 


1.1 


42. 


.014543 


-3. 


.5284089 


1.01 




1.2 


42. 


.014346 


-3 


.5287867 






1.3 


42. 


.009577 


-3 


.5304728 






2.1 


42. 


.012115 


-3 


.5330117 


1.01 




2.2 


42. 


.011774 


-3 


.5327867 






2.3 


42. 


.010140 


-3 


.5313950 






3.1 


42. 


.009979 


-3 


.5331978 


2.55 




3.2 


42. 


.009479 


-3 


.5331172 






3.3 


42. 


.012463 


-3 


.5352256 






4.1 


42. 


.011771 


-3 


.5352367 


2.55 




4.2 


42. 


.009213 


-3 


.5339200 






4.3 


42. 


.009082 


-3 


.5334061 






5.1 


42. 


.013618 


-3 


.5263012 


6.03 




5.2 


42. 


.019169 


-3 


.5328972 






6.1 


42. 


.017633 


-3 


.5313868 




1.93±0.15 


6.2 


42. 


.013918 


-3 


.5331647 






6.3 


42. 


.008825 


-3 


.5280806 






6.4 


42. 


.015347 


-3 


.5266792 






7.1 


42. 


.016975 


-3 


.5238458 




>4.0 


7.2 


42. 


.014894 


-3 


.5231097 






7.3 


42. 


.013315 


-3 


.5229153 






8.1 


42. 


.015123 


-3 


.5234406 




2.19±0.29 


8.2 


42. 


.014163 


-3 


.5232253 






9.1 


42. 


.016088 


-3 


.5336001 




>4.0 


9.2 


42. 


.016702 


-3 


.5331010 






9.3 


42. 


.016036 


-3 


.5264248 






9.4 


42. 


.007825 


-3 


.5278831 







Table 1. Multiply imaged systems considered in this work. 
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ment with previous works, we find that Abell 383 is well 
described by an unimodal mass distribution which presents 
a small eccentricity. This indicates a well relaxed dynamical 
state. We note that the galaxy scale perturbers all present 
a scale radius which is smaller than the scale radius inferred 
for isolated field galaxies, in agreement with the tidal strip- 
ping scenario. 

The Lenstool software does explore the parameter 
space using a Monte Carlo Markov Chain sampler. At the 
end of the optimization, we have access to these MCMC re- 
alizations from which we can draw statistics and estimate 
error bars. For each realization, we build a two dimensional 
mass map. All these mass maps are then used to compute the 
mean mass map and the corresponding covariance. These in- 
formation will then be used in the joint fit. 



2.4 The case of system 6 

Following IZitrin et all l |2011bl ), we have included the four 
images constituting system 6 in our analysis and our model 
reproduces it correctly. However, image 6.2 exhibits a seri- 
ous problem of symmetry. This is the first time we observe 
such a problem in a strong lensing configuration and we do 
not have any clear explanation for this. As can be seen on 
Fig. [2] images 6.1 and 6.3 are constituted by a bright spot 
linked to a faint tail. The symmetry between these images 
is correct. Image 6.2 does exhibit two tails. One may be 
tempted to naturally associate Tail 1 with image 6.2. How- 
ever, lensing symmetry requires the associated tail to point 
to the west. Tail 2 could be the associated feature, and an 
eventual perturber may bend this tail a bit north. 

Another explanation could be that a undetected dark 
matter substructure does contract the expected tail into im- 
age 6.2. 



3 X-RAY DATASETS AND ANALYSIS 

The cluster Abell 383 is a luminous cluster at redshift z — 
0.189, which exhibits several indications of a well relaxed 
dynamical state, for instance the absence of evident sub- 
structures and a central X-ray surface brightness peak, as- 
sociated with a cool core. The global (cooling-core corrected) 
temperature T ew has been estimated to be T ow = 4.17±0.10 
keV and an abundance of 0.43 ± 0.06 solar value ( t|3.1|l . We 
classify this cluster as a strong cooling core source (SCC, 
iMorandi fc Ettoril feoOT). i.e. the central cooling time £ C ooi 
is less than the age of the universe £ agc ,z at the cluster red- 
shift (t C ooi/iagc,z < 0.1): we estimated a t coo i — 7 X 10 8 yr. 
As other SCC sources, Abell 383 shows a very low central 
temperature (~ 2.7 keV) and a strong spike of luminosity in 
the brightness profile. The temperature profile is very regu- 
lar suggesting a relaxed dynamical state (see upper panel of 
Fig. H|). 

Des cription of the X-ray analys is methodology can be 
found in IMorandi et all (|2007l . l2010f ). Here we briefly sum- 
marize the most relevant aspects of our data reduction and 
analysis of Abell 383. 



3.1 X-ray data reduction 

We performed our X-ray analysis on three datasets retrieved 
from the NASA HEASARC archive (observation ID 2320, 
2321 and 524) with a total exposure time of approximately 
50 ks. We summarize here the most relevant aspects of the 
X-ray data reduction procedure for Abell 383. Two obser- 
vations have been carried out using ACIS-I CCD imaging 
spectrometer (ID 2320 and 524), one (ID 2321) using ACIS- 
S CCD imaging spectrometer. We reduced these observa- 
tions using the CIAO software (version 4.3) distributed by 
the Chandra X-ray Observatory Center, by considering the 
gain file provided within CALDB (version 4.4.3) for both 
the observations telemetered in Very Faint mode (ID 2320 
and 524) and Faint mode (ID 2321). 

We reprocessed the level- 1 event files to include the 
appropriate gain maps and calibration products. We used 
the acis_process_events tool to check for the presence of 
cosmic-ray background events, correct for spatial gain varia- 
tions due to charge transfer inefficiency and re-compute the 
event grades. Then we have filtered the data to include the 
standard events grades 0, 2, 3, 4 and 6 only, and therefore we 
have filtered for the Good Time Intervals (GTIs) supplied, 
which are contained in the f ltl .fits file. We then used the 
tool dmextract to create the light curve of the background. 
In order to clean the datasets of periods of anomalous back- 
ground rates, we used the deflare script, so as to filter 
out the times where the background count rate exceed ±3<r 
about the mean value. Finally, we filtered ACIS event files 
on energy selecting the range 0.3-12 keV and on CCDs, so 
as to obtain an level-2 event file. 



3.2 X-ray spatial and spectral analysis 

We measure the gas density profile from the surface bright- 
ness recovered by a spatial analysis, and we infer the pro- 
jected temperature profile by analyzing the spectral data. 

The X-ray images were extracted from the level-2 event 
files in the energy range (0.5 — 5.0 keV), corrected by the 
exposure map to remove the vignetting effects, masking out 
point sources and by rebinning them of a factor of 4 (1 
pixel= 1.968 arcsec). 

We determined the centroid (x c ,y c ) of the surface 
brightness by locating the position where the X and Y 
derivatives go to zero, which is usually a more robust de- 
termination than a center of mass or fitting a 2D Gaussian 
if the wings in one direction are affected by the presence of 
neighboring substructures. We checked that the centroid of 
the surface brightness is consistent with the center of the 
BCG, the shift between them being evaluated to be ~ 1.5 
arcsec: note that the uncertainty on the X-ray centroid es- 
timate is comparable to the applied rebinning scale on the 
X-ray images. 

The three X-ray images were analyzed individually, in 
order to calculate the surface brightness for each dataset. 

The spectral analysis was performed by extracting the 
source spectra in circular annuli of radius r* n around the 
brightest cluster galaxy (BCG). We have selected n* = 7 
annuli out to a maximum distance i? S pcc = 1056 kpc in order 
to have a number of net counts of photons of at least 2000 
per annulus. All the point sources have been masked out 
by both visual inspection and the tool celldetect, which 
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Clump 


R.A. 


Decl. 


e 


e 


a (kms 1 ) 


r (arcsec) 


Halo 
cD 


0.98±0.29 

[o.o] 


0.54±0.45 
[0.0] 


0.18±0.02 
[0.105] 


112.0±1.3 
[101] 


907±20 
332±« 


15.6±1.3 
>15 


PI 


[14.7] 


[-16.9] 


[0.185] 


[157] 


179±f 


o 7 +10.3 
°-'-0.6 


P2 


[9.0] 


[-2-0] 


[0.589] 


[184] 


<150 


<6.0 


elliptical galaxy 










<165 


<4.0 



Table 2. Mass model parameters. Coordinates are given in arcseconds with respect to the cD Galaxy. The ellipticity e is the one of the 
mass distribution, expressed as (a 2 — b 2 )/ (a 2 + b 2 ) . Error bars correspond to la confidence level as inferred from the MCMC optimization. 
When the posterior probability distribution is not Gaussian, we report the mode and asymmetric error bars. Values into brackets are not 
optimized. The radius reported here corresponds to the core radius for the halo and to the scale radius for the other mass components. 




provide candidate point sources. Then we have calculated 
the redistribution matrix files (RMF) and the ancillary re- 
sponse files (ARF) for each annulus: in particular we used 
the CI AO specextract tool both to extract the source and 
background spectra and to construct the ARF and RMF of 
the annuli. 

For each of the n* ann uli the spectra have been ana- 
lyzed by using the XSPEC |A rnaudl ll996) packag e, by si- 
multaneously fitting an a bsorbed MEKAL model (|Kaastral 
ll992l ; lLiedahl et al.l ll995) to the three observations. The fit 
is performed in the energy range 0.6-7 keV by fixing the 
redshift at z = 0.189, and the photoelectric absorption at 
the galactic value. For each of the n* annuli the spectra 
we grouped the photons into bins of 20 counts per energy 
channel and applying the x 2 -statistics. Thus, for each of the 
annuli, the free parameters in the spectral analysis were the 
normalization of the thermal spectrum K\ oc J nj? dV, the 
emission- weighted temperature T* ro j ; , and the metallicity 
Zi. 

The three observations were at first analyzed individ- 
ually, to assess the consistency of the datasets and to ex- 
clude any systematic effects that could influence the com- 
bined analysis. We then proceeded with the joint spectral 
analysis of the three datasets. 

The background spectra have been extracted from re- 
gions of the same exposure for the ACIS-I observations, for 



which we always have some areas free from source emission. 
We also checked for systematic errors due to possible source 
contamination of the background regions. Conversely, for the 
ACfS-S observation we have considered the ACIS-S3 chip 
only and we used the ACIS "blank-sky" background files. 
We have extracted the blank sky spectra from the blank- 
field background data sets provided by the ACIS calibration 
team in the same chip regions as the observed cluster spec- 
tra. The blank-sky observations underwent a reduction pro- 
cedure comparable to the one applied to the cluster data, 
after being reprojected onto the sky according to the obser- 
vation aspect information by using the reproject_events 
tool. We then scaled the blank sky spectrum level to the 
corresponding observational spectrum in the 9-12 keV inter- 
val, where very little cluster emission is expected. We also 
verified that on the ACIS-I observations the two methods of 
background subtraction provide very similar results on the 
fitting parameters (e.g. the temperature). 



4 THREE-DIMENSIONAL STRUCTURE OF 
GALAXY CLUSTERS 

The lensing and the X-ray emission both depend on the 
properties of the DM gravitational potential well, the for- 
mer being a direct probe of the two-dimensional mass map 
via the lensing equation and the latter an indirect proxy of 
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the three-dimensional mass profile through the HE equation 
applied to the gas temperature and density. In order to infer 
the model parameters of both the IC gas and of the underly- 
ing DM density profile, we perform a joint analysis of SL and 
X-ray data. We briefly outline the methodology in order to 
infer physical properties in triaxial galaxy clusters: (1) We 
start with a generalized Navarro, Frenk and White (gNFW ) 
triaxial model of the DM as described in ljing fe Sutol 12002). 
which is representative of the total underlying mass distri- 
bution and depends on a few parameters to be determined, 
namely the concentration parameter C, the scale radius i? s , 
the inner slope of the DM 7 , the two axis ratios (r]DM,a 
and r/QM.b) and t h e Eul er angles tp, 9 and <j> (2) following 
iLee fc Sutol l|2003l , | 20041 ). we recover the gravitational po- 
tential (Equation [8| and two-dimensional surface mass S 
(Equation 1 14[) of a dark halo with such triaxial density pro- 
file; (3) we solve the generalized HE equation, i.e. including 
the non-thermal pressure P n t (Equation 1 1 Of) . for the den- 
sity of the IC gas sitting in the gravitational potential well 
previously calculated, in order to infer the theoretical three- 
dimensional temperature profile T gas in a non-parametric 
way; (4) we calculate the surface brightness map Sx related 
to the triaxial ICM halo (Equation I13[) ; and (5) the joint 
comparison of T gas with the observed temperature, of Sx 
with observed brightness image, and of £ with the observed 
two-dimensional mass map gives us the parameters of the 
triaxial ICM and DM density model. 

Here we briefly summarize the major findings of 
iMorandi et al.l (|2010h for the joint X-ray+Lensing analysis 
in order to infer triaxial physical properties, as well as the 
improvements added in the current analysis; additional de- 
tails can be found in IMorandi et all (|2007l , |2010| . l2011al ibh. 

We start by describing in Sect. l4~Tl the adopted geom- 
etry, the triaxial DM density and gravitational potential 
model, focusing on the relation between elongation of ICM 
and DM ellipsoids. We describe the relevant X-ray and lens- 
ing equations in £14.21 then we discuss how to jointly combine 
X-ray and lensing data in H4.3I 



4.1 ICM and DM triaxial halos 

Extending our previous works, in the present study we allow 
the DM and ICM ellipsoids to be orientated in a arbitrary 
direction on the sky. We introduce two Cartesian coordinate 
systems, x = (x,y,z) and x' = (x',y',z'), which represent 
respectively the principal coordinate system of the triaxial 
dark halo and the observer's coordinate system, with the 
origins set at the center of the halo. We assume that the z'- 
axis lies along the line of sight direction of the observer and 
that the x' , y' axes identify the direction of West and North, 
respectively, on the plane of the sky. We also assumed that 
the x, y, z-axes lie along the minor, intermediate and major 
axis, respectively, of the DM halo. We define tp, 9 and <f> as 
the rotation angles about the x, y and z axis, respectively 
(see Figure[3]). Then the relation between the two coordinate 
systems can be expressed in terms of the rotation matrix M 
as 



x' = Mx, 



(1) 




Figure 3. The orientations of the coordinate systems. The 
Cartesian axes (x,y,z) represent the DM and ICM halo prin- 
cipal coordinate system while the axes (x',y',z') represent the 
observer's coordinate system with z'-axis aligned with the line 
of sight (l.o.s.) direction. We define tp, 8 and <j> as the rotation 
angles about the x, y and z axis, respectively. The labels N and 
W indicate the position of the North and West, respectively, on 
the plane of the sky. 



M x (tp),M y (9),M z (cp) with Euler angles ip,9,<p, and it is 
given by: 



M = M x (ip)#M y (0)#M z (<p) . 



(2) 



where 



My (6) 



M z (<f>) = 



1 

cos tp — sin tp 
sin tp cos tp 

cos 9 sin 9 

1 
— sin 9 cos 9 



(3) 



COS ( 

sin <; 




— sin <p 
cos (p 
1 



Figure [3] represents the relative orientation between the 
observer's coordinate system and the halo principal coordi- 
nate system. 

In order to parameterize the cluster mass distribution, 
we consider a triaxial generalized Navarro, Frenk & White 
model gNFW (|jing fc Sutoll2002D : 



p(R) 



Sopc,z 



{R/RsT (1 + R/Rs) 



3-7 ' 



(4) 



where M represents the orthogonal matrix corre- 
sponding to counter-clockwise/right-handed rotations 



where R s is the scale radius, 8c is the dimensionless char- 
acteristic density contrast with respect to the critical den- 
sity of the Universe pc,z at the redshift z of the cluster, 



© 0000 RAS, MNRAS 000, 000-000 



Triaxiality and non-thermal pressure in A3 83 7 



and 7 represents the inner slope of the density profile; 
pc,z = 3H (z) 2 /8ivG is the critical density of the universe 
at redshift z, H z = E z Ho, E z = [« M (1 + zf + Ov] 1/2 , and 

200 Cf (5) 



5c = 



3 F(C, 7 ) ' 



where C = -R 200 / Rs is the concen tration parameter. F(C, 7) 
is defined as i|Wvithe et al.ll200ll ): 



F(C, 7 ) = [° s 2 -<(i + sy- 3 ds. 

Jo 



(6) 



The radius R can be regarded as the major axis length 
of the iso-density surfaces: 



2 / x 2 v 2 z 2 
R = c I — + ^ + 



(a ^ b ^ c). 



(7) 



We also define rjDM,a — a/c and T)nM,b — b/c as the minor- 
major and intermediate-major axis ratios of the DM halo, 
respectively. 

The gravitational potential of a dark halo with the tri- 
axial density profile ( Equation [4]) can be writte n as a com- 
plex implicit integrals dBinnev fc T rcmainc 1987). While nu- 
merical integration is requi red in genera l to ob tain the triax- 
ial gravitational potential, iLee fe Sutol (|2003l ') retrieved the 
following approximation for the gravitational potential $ 
under the assumption of triaxial gNFW model for the DM 
(Equation Q}: 



*(u) ~C Fi(v) + Co ^±^F 2 (u) 



+C 



> + e c cos 



(8) 



-F 3 (u), 



with u = r/Rs, Co = 47rG<5 c p C ritPs , and the t hree functions, 
Fi(u) , F2 (it), and ^3(14) has been defined in iMorandi et all 
(2010), eb (eb) and e c (e c ) are the eccentricity of DM (IC 
gas) with respect to the major axis ( e.g. ej, = \J\ — (b/c) 2 ). 

The work of ILee fc Sutol |2003l ) showed that the iso- 
potential surfaces of the triaxial dark halo are well ap- 
proximated by a sequence of concentric triaxial distribu- 
tions of radius Pi cm with different eccentricity ratio. For 
Pi cm it holds a similar definition as R (Equation 0), but 
with eccentricities e;, and e c . Note that eb = tb(fib, u, 7) and 
e e = e c (e c ,u, 7), unlike the constant et, e c for the adopted 
DM halo profile. In the whole range of u, ebl&b (cc/e c ) is less 
than unity (~ 0.7 at the center), i.e., the intracluster gas is 
altogether mor e sphe rical than the underlying DM halo (see 
IMorandi etafl (|2010h for further details). 

The iso-potential surfaces of the triaxial dark halo co- 
incide also with the iso-density (pressure, temperature) sur- 
faces of the intracluster gas. This is si mply a direct con- 
seque nce of the X-ray shape theorem (|Buote fc Canizaresl 
1 19941 ); the HE equation (JlOj) yields 



VP x V$ = VjOgas x V$ = 0. 



4.2 X-ray and lensing equations 



(9) 



For the X-ra y analysis we rely on a generalization of the 
HE equation l|Morandi et al.ll2011bh . which accounts for the 
non-thermal pressure P n t and reads: 



where p gas is the gas mass density, $ is the gravitational po- 
tential, Ptot = Pth + P n t, and with the non-thermal pressure 
of the gas P n t assumed to be constant fraction £ of the total 
pressure Ptot, i.e. 



Pnt = £Ptc 



(11) 



Note that X-ray data probe only the thermal component of 
the gas P t h = riekPgas, k being the Boltzmann constant. 
From Equations (|10|) and (|11|) we point out that neglecting 
P n t (i.e. Ptot = Pth) systematically biases low the determi- 
nation of cluster mass profiles roughly of a factor £. 

To model the density profile in the triaxial ICM halo, 
we use the following fitting function: 

MPicm) = no (R icm /r c )- S (l + Rf crn /r 2 c )- 3/2s+s/2 (12) 

with parameters (no, r c , e, 8). We computed the theoretical 
tree-dimensional temperature T gas by numerically integrat- 
ing the equation of the HE (Equation [TDJ, assuming triaxial 
geometry and a functional form of the gas density given by 
Equation (JT2J . 

The observed X-Ray surface brightness Sx is given by: 



Sx = 



47r(l + z) 



n e n p dz 



(13) 



where A(T* ro j,Z) is the cooling function. Since the pro- 
jection on the sky of the plasma emissivity gives the X- 
ray surface brightness, the latter can be geometrically fit- 
ted with the model n e (Pi cm ) of the assumed distribution 
of the gas density (Equation I12|) by applying Equation (|13p . 
This has been accomplished via fake Chandra spectra, where 
the current model is folded through response curves (ARF 
and RMF) and then added to a background file, and with 
absorption, temperature and metallicity measured in that 
neighboring ring in the spectral analysis (i )3.2|l . In order to 
calculate A(Tp rn; , Z) , we adopted a MEKAL model (|Kaastral 
Il992l; iLiedahl et al]|l995l ) for the emissivity. 

For the lensing analysis the two-dimensional surface 
mass density S can be expressed as: 



p{R)dz' 



(14) 



VPtot = -pga S V$ 



(10) 



We also calculated the covariance matrix C among all 
the pix els of the observed surface mass (see IMorandi et aD 
(|2011bT ) for further details). 



4.3 Joint X-ray+lensing analysis 

The lensing and the X-ray emission both depends on the 
properties of the DM gravitational potential well, the former 
being a direct probe of the projected mass profile and the 
latter an indirect proxy of the mass profile through the HE 
equation applied on the gas temperature and density. In this 
sense, in order to infer the model parameters, we construct 
the likelihood performing a joint analysis for SL and X-ray 
data, to constrain the properties of the model parameters of 
both the ICM and of the underlying DM density profile. 

The system of equations we simultaneously rely on in 
our joint X-ray+Lensing analysis is: 
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Tgas = T gas (C, Rs, 7, T]DM,a, T]DM,b, i> , 9, <f>, no, T c , S, S, f) 

Sx = Sx(C, R B ,j,r/DM,a,VDM,b,ip,6,4>,no,r c ,£,5) (15) 

where the three-dimensional model temperature T gas is re- 
covered by solving equation (|10[) and constrained by the ob- 
served temperature profile, the surface brightness is recov- 
ered via projection of the gas density model (Equation 1 13[) 
and constrained by the observed brightness, and the model 
two-dimensional mass £ is recovered via Equation (| 14[) and 
constrained by the observed surface mass. 

The method works by constructing a joint X- 
ray+Lensing likelihood C oc exp(— x 2 /2): 

X = Xx,t + Xx,s + xL s (16) 

with Xx,t, xl,s and xLs being the y 2 coming from the X- 
ray temperature, X-ray brightness and lensing data, respec- 
tively. 

For the spectral analysis, Xx,t ls equal to: 

2 _ \ ~~ (?proj,i " ?proj,i) 

Xx,T - 2^ Z2 ( n > 

i-1 T" ■ ■ 

Proja 

^proj.i being the observed projected temperature profile in 
the ith circular rin g and r pr oi.i the azim uthally-averaged 
projection (following Ma zzotta et al.ll2004h of the theoretical 
three-dimensional temperature Tgas; the latter is the result 
of solving the HE equation, with the gas density n e (i?i cm ). 
For the X-ray brightness, x x ,s reads: 

2 \r^sr^ (Sx,i - S x ,i) 2 , . 

xx,s = 2^ -2 ( 18 ) 

j i=l a s,i 

with Sx,i and S x ,i theoretical and observed counts in the ith 
pixel of the jth image. Given that the number of counts in 
each bin might be small (< 5), then we cannot assume that 
the Poisson distribution from which the counts are sampled 
has a nearly Gaussian shape. The standard deviation (i.e., 
the square-root of the variance ) for this low-count case has 
been derived bv lGehrelsl (|l986D : 

a s , t = 1 + Js^t + 0.75 (19) 

which has been proved to be accurate to approximately one 
percent. Note that we added background to Sx,i as mea- 
sured locally in the brightness images, and that the vi- 
gnetting has been removed in the observed brightness im- 
ages. 

For the lensing constraints, the lensing contribution is 
X?ens = [S-E*] t C- 1 [E-E*] , (20) 

where C is the covariance matrix of the 2D projected mass 
profile from strong lensing data, S* = (Si, SJ, EJj*) are 
the observed measurements of the projected mass map in the 
ith pixel, and £ is the theoretical 2D projected mass within 
our triaxial DM model. Note that we removed the central 25 
kpc of the 2D projected mass in the joint analysis, to avoid 
the contamination from the cD galaxy mass. 

The probability distribution function of model param- 
eters has been evaluated via Markov Chain Monte Carlo 
(MCMC) algorithm. Errors on the individual parameters 



have been evaluated by considering average value and stan- 
dard deviation on the marginal probability distributions of 
the same parameters. 

So we can determine the physical parameter of 
the cluster, for example the 3D temperature T gas = 
T gas (C, R B ,y,riDM,a,VDM,b,4>,8,(f),no,r c ,£,6,£), the shape 
of DM and ICM, just by relying on the HE equation and 
on the robust results of the hydrodynamical simulations of 
the DM profiles. In Fig. [3] we present an example of a joint 
analysis for T gas , Sx and S: for Sx and £ the ID profile has 
been presented only for visualization purpose, the fit being 
applied on the 2D X-ray brightness/surface mass data. Note 
that in the joint analysis both X-ray and lensing data are 
well fitted by our model, with a Xrcd — 1-23(65837 degrees 
of freedom). 



5 RESULTS AND DISCUSSION 

In the previous section we showed how we can determine 
the physical parameters of the cluster by fitting the avail- 
able data based on the HE equation and on a DM model 
that is based on robust results of hydrodynamical cluster 
simulations. In this section we present our results and dis- 
cuss their main implications. We particularly focus on the 
implications of our analysis for the determination of the full 
triaxiality, viability of the CDM scenario, the discrepancy 
between X-ray and lensing masses in Abell 383, and the 
presence of non-thermal pressure. 

5.1 Best-fit parameters 

In table [3] we present the best-fit model parameters for our 
analysis of Abell 383. Our work shows that Abell 383 is a 
triaxial galaxy cluster with DM halo axial ratios r)r>M,a = 

0. 55±0.06 and r]r>M,b = 0.71±0.10, and with the major axis 
inclined with respect to the line of sight of 9 — 21.1 ± 10.1 
deg. Note that these elongations are statistically significant, 

1. e. it is possible to disprove the spherical geometry assump- 
tion. We also calculated the value of M200 ■ 

800?r n 3 
M200 = —jr—riDM,ar)DM,bft<iaoPc,z , (21) 

We have: M 200 = 7.3 ± 0.5 x 10 14 A/ Q and R 200 = 2401.1 ± 
166.8 kpc. 

The axial ratio of the gas is ??gas, a ~ 0.70—0.82 and 
^gas,b ~ 0.80—0.88, moving from the center toward the X- 
ray boundary. 

Another important result of our work is the need for 
non-thermal pressure support, at a level of ~10%. 

We also observe that the small value of the inclination 
of the major axis with respe ct to the line of sight is in agree- 
ment with the predictions of lOguri fc Blandfordl {2009), who 
showed that SL clusters with the largest Einstein radii con- 
stitute a highly biased population with major axes prefer- 
entially aligned with the line of sight thus increasing the 
magnitude of lensing. 

In fig [5] we present the joint probability distribution 
among different parameters in our triaxial model. For ex- 
ample we observe that there is a positive (negative) linear 
correlation between 7 — 9(C). This proves that e.g. the incli- 
nation with respect to the line of sight, ad more in general 



© 0000 RAS, MNRAS 000, 000-000 



Triaxiality and non-thermal pressure in A3 83 9 




200 



400 600 
R (kpc) 



800 



1000 



Figure 4. Example of the joint analysis for T gas , Sx and S. In the upper panel we display the two quantities which enter in the X-ray 
analysis (Equation |17> : the observed spectral projected temperature T*- m (big points with errorbars) and the theoretical projected 
temperature T pro j m (diamonds). We also show the theoretical 3D temperature (points), which generates T pro j iIn through convenient 
projection techniques. In the middle panel we display the two quantities which enter in the X-ray brightness analysis (Equation llSII : the 
observed surface brightness profile S x (points with errorbars) and the theoretical one Sx (solid line). In the lower panel we display the 
two quantities which enter in the lensing analysis (Equation H20H ): the observed surface mass profile S* (points with error bars) and the 
theoretical one S (solid line). Note that for surface brightness (surface mass) the ID profile has been presented only for visualization 
purpose, the fit being applied on the 2D X-ray brightness (lensing) data. Moreover, for the surface brightness we plotted data referring 
to the observation ID 2320. 



Table 3. Best-fit model parameters of Abell 383. 



c 


4.76 ±0.51 


Rs (kpc) 


511.2 ±73.6 


7 


1.02 ±0.06 


VDM,a 


0.55 ±0.06 


VDM,b 


0.71 ±0.10 


tp (deg) 


-13.6 ±5.5 


6 (deg) 


21.1 ± 10.1 


<t> (deg) 


-16.9 ± 15.9 


n (cm -3 ) 


0.063 ± 0.003 


r c (kpc) 


26.4 ± 1.7 


£ 


0.55 ±0.01 


8 


0.02 ±0.01 


t 


0.11 ±0.05 



the geometry, is important in order to characterize the in- 
ner slope of the DM (and the concentration parameter), and 
therefore to assess the viability of the standard cosmological 
framework. 



We also compared the azimuthal angle (f> — 103.8 ± 9.0 
deg and the eccentricity on the plane of the sky (e = 0.08 ± 
0.03), with the values on the total two-dimensional mass 
from the analysis from Lenstool (0 = 111 ± 2 deg and 
e = 0.05 ± 0.01): note the good agreement. For t he me thod 
to recover e and d> we remand to iMorandi etlll (|2010T ). 



5.2 Probing the inner DM slope and implications 
for the viability of the CDM scenario 

A central prediction arising from CDM simulations is that 
the density profile of DM halos is universal across a wide 
ran ge of mass scales from dwarf galaxies to clusters of galax- 
ies (|Navarro et all 1 19971 ): within a scale radius, r B , the DM 
density asymptotes to a shallow powerlaw trend, Pdm(t) oc 
r -7 , with 7=1, while external to r s , Pdm(t) oc r~ 3 . 
Nevertheless, the value of the logarithmic inner slope o is 
still a matter of debate. Several s tudie s suggested steeper 
central cusp dMoore et al.l 1 19981 . 1 19991): with one recent 
suite of high-resolution sim ulations iDiemand et"aL I J2004D 
argued for 7 = 1.16 ± 014. i Tavlor fc Navarro! (|200ll ) and 
iDehnen fc McLaughlin! (|2005F T developed a general frame- 
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work where there is the possibility that 7 might be shallower 
than 1 at small radii, while other studies suggested th e lack 
of cusp in the innermost regions (jNavarro et al.l [20041 ). 

Conventional CDM simulations only include collision- 
less DM particles, neglecting the interplay with baryons. It 
is not clear how the inclusion of baryonic matter would af- 
fect the DM distribution: although their impact is thought 
to be smaller in galaxy clusters than in galaxies, baryons 
may be crucially important on scales comparable to the ex- 
tent of typical brightest cluster galaxies and may thus alter 
the dark matter distribution, particularly the inne r slope 
( San d et al 1 120081 : Is ommer-Larsen fc Lim ousin 200||). With 
this respect, the cooling of the gas is expected to steepen the 
DM distribution via adiabatic compression (jGnedin et al.l 
l2004h . Dynamical friction between infalling baryons and the 
host halo may counter and surpass the effects of adiabatic 
contraction in galaxy clusters, by driving energy and angular 
momentum out of the cluster core, thu s softening an origi - 
nally cuspy profile (|E1-Zant et al.ll2004h . [Mead et all (|2010j ) 
shows that AGN feedback suppresses the build-up of stars 
and gas in the central regions of clusters, ensuring that the 
inner density slope is shallower and the concentration pa- 
rameter smaller than in cooling-star formation haloes, where 
C and 7 are similar to the dark-matter only case. 

An observational verification of the CDM predictions, 
via a convincing measurement of 7 over various mass scales, 
has proved controversial despite the motivation that it of- 
fers a powerful test of the CDM paradigm. Observation- 
ally, efforts have been put on probing the central slope 
7 of the underlying dark matter distribution , through X- 
ray (jEttori et all 120021 : lArabadiis et al.l 120021 : iLewis et al.l 



20031: IZappacosta et all 120061^1: l e nsing dlYson et al 



Smith et al.l 


200 ll: 


Gavazzi et al. 


120031: 



Dahle et al. 120031; Sand et all 



l2003l: Gavazzi 2005: ISand eTa! 



iLimousin et al] 12008 



200. 




or dynamics 

ijKelson et al.l2002i : lBiviano fc Salu cci 2006). However, anal- 
yses of the various data sets have given conflicting values of 
7, with large scatter from one cluster to another, and many 
of the assumptions used have been questioned. Indeed these 
determinations rely e.g. on the standard spherical modeling 
of galaxy clusters: possible elongation/flattening of the clus- 
ters along the line of sight has been proved to affect the esti- 
mated values of 7 (jMorandi et al.ll20ld . l2011ah . Another ma- 
jor observational hurdle is the importan ce of separating th e 
baryonic and non-baryonic components (|Sand et al]|200a ). 

With this perspective, one of the main result of the 
presented work is to measure a central slope of the DM 
1.02 ± 0.06 by accounting explicitly for the 3D structure 
for Abell 383: this value is in agree ment with the CDM pre- 
dictions from lNavarro et al] (|l997f ) (i.e. 7 = 1). Yet, a com- 
parison with such simulations, which do not account for the 
baryonic component, merits some caution. Note that in this 
analysis we did not subtract the mass contribution from the 
X-ray gas to the total mass. Nevertheless the contribution of 
the gas to the total matter is small: the measured gas frac- 
tion is 0.03 — 0.05 in the spatial range 25 — 400 kpc, and we 
checked that the slope of the density profile is very similar 
to that of the DM beyond a characteristic scale r c ~ 20 — 30 
kpc, a se lf-similar property of the gas common to all the SCC 
sources (jMorandi fc Ettorill2007l ) . This suggests that our as- 
sumption to model the total mass as a gNFW is reliable, 
and therefore a comparison with DM-only studies is ten- 



able. S imil ar conclusions have been reached bylBradac et al] 
(|2008h and lSommer-Larsen fc Limousin! ij2009l ). 

It is also possible that the DM particle is self- 
interacting: this would lead to shallower DM density 
profiles with respect to the standard CDM scenario 
ijSpergel fc Steinhardtl |2000| ). Nevertheless, the lack of a 
flat core allows us to put a conservative upper limit 
on the dark matter p a rticle scattering cross section: in- 
deed e.g. lYoshida et al] l|2000l ) simulated cluster-sized halos 
and found that relatively small dark matter cross-sections 
(crdm = 0.1cm 2 g" 1 ) are ruled out, producing a relatively 
large (40/t -1 kpc) cluster core, which is not observed in our 
case study. 

iMorandi et al] (|2010l . l2011al ) analyzed the galaxy clus- 
ters Abell 1689 and MACS J1423. 8+2404 in our triaxial 
framework: we found that 7 is close to the CDM predictions 
(7 = 1) for these clusters, in agreement with the results in 
the present pa per. 

Recently, iNewman et al] (|201ll ) combined strong and 
weak lensing constraints on the total mass distribution with 
the radial (circularized) profile of the stellar velocity disper- 
sion of the cD galaxy. They adopted a triaxial DM distri- 
bution and axisymmetric dynamical models, and they also 
added hydrostatic mass estimates from X-ray data under the 
assumption of spherical geometry as a further constraint. 
The proved that the logarithmic slope of the DM density at 
small radii is 7 < 1.0 (95% confidence level), in disagree- 
ment with the present work. Note that, although the geo- 
m etrical model employe d in the work is more accurate than 
in lNewman et afl (|201ll ) (see H5.3l for a discussion about the 
possible systematics involved in their analysis), we do not 
have constraints from stellar velocity dispersion of the cD 
galaxy: therefore, a strict comparison between the two works 
merits some caution. 

With this perspective, we plan to extend our work in 
the future by including in our triaxial framework the stellar 
velocity dispersion of the cD galaxy, which is essential for 
separating the baryonic mass distribution in the cluster core: 
this will allow a more tight determination of the inner slope 
of the DM on scales down to ~ a few kpc, while we removed 
the central 25 kpc of the 2D projected mass, to avoid the 
contamination from the cD galaxy, though we checked that 
there is a very little dependence of the physical parameters 
on the radius of the masked region. 

We also plan to collect data for a larger sample of clus- 
ters in order to characterize the physical properties of the 
DM distribution. Indeed we proved that the inner slope of 
the DM is in agreement with the theoretical predictions just 
for a few clusters, while a broader distribution of inner slopes 
might exist in nature than that predicted by pure dark mat- 
ter simulations, pote ntially depending e.g. on the merge r 
history of the cluster (jNavarro et alj|2004l ; lGao et al.ll201ll ). 

The value of the concentration parameter C = 4.76 ± 
0.51 is in agreement w ith the theoretica l expectation from 
N-body simulations of iNeto et al] (|2007l ). where C ~ 4 at 
the redshift and for the virial mass of Abell 3 83, and with 
an intr insic scatter of ~ 20 per cent. Recently. IZitrin et al] 
lj2011bl ) reported a high value of the concentration CW = 
8.591J5 20^0 41 Yl& a joint strong and weak lensing analysis 
under the assumption of spherical geometry. Their concen- 
tration parameter refers to the virial radius (we refer to 
R20Q) and includes the statistical followed by the system- 
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atic uncertainty. They argue that Abell 383 lies above the 
standard C — M relation, putting some tension with the pre- 
dictions of the sta ndard model. Our de termination of C is in 
disagreement with lZitrin et al] l|2011bl ). as for other clusters 
with prominen t strong l ensing features (e.g. Abell 1689, see 
iMorandi et al.l (|2011al w), due to the improved joint X-ray 



and strong l ensing triaxial analysis we implemented. Yet, 
IZitrin et al.l l|2011bT ) observed that, though even with effects 
of projections taken into account in numerical simulations, 
there seems to be some marginal discrepancy from ACDM 
predictions of simulated adiabatic clusters extracted from 
the MareNostrum Univ erse in terms of lensing efficiency 
l|Meneghetti et all |201 lh . This departure might be due to 
selection biases in t he observed subsample at this redshift 
l|Horesh et alj|201ll ) or to physical effects that arise at low 
redshift and enhance the lensing efficiency and not included 
the previous simul ations, as the effect of b aryons on the DM 
profile in clusters (jBarkana fc Loebll201Cl 1. 



5.3 Implications on X-ray/Lensing mass 
discrepancy 

Here we outline our findings in probing the 3D shape of 1CM 
and DM, and the systematics involved in using a standard 
spherical modeling of galaxy clusters. In particular, we will 
focus on the long-standing discrepancy between X-ray and 
lensing masses on clusters, showing that this is dispelled if 
we account explicitly for a triaxial geometry. 

There are discrepancies between cluster masses deter- 
mined based on gravitational lensing and X-ray observa- 
tions, the former being significantly higher than the lat- 
ter in clusters with pro minent lensing features. Indeed 
lOeuri fc Blandfordl [|2009h showed that SL clusters with the 
largest Einstein radii constitute a highly biased population 
with major axes preferentially aligned with the line of sight 
thus increasing the magnitude of lensing. Given that lensing 
depends on the integrated mass along the line of sight, ei- 
ther fortuitous alignments with mass concentrations that are 
not physically related to the galaxy cluster or departures of 
the DM halo from spherical symmetry can bias up t he three- 
dime nsional with mass respect to the X-ray masses (|Gavazzil 
2005); on the other hand, X-ray-only masses hinges on the 
goodness of the strict HE assumption: the presence of bulk 
motions in the gas can bias l ow the three-dimensiona l mass 
profile between 5% and 20% (jMeneehetti et alJl201oh . 

In Figure [6] we compare the 3D mass enclosed within a 
spherical apertures of radius R for lensing, X-ray-only data 
and from a joint X-ray+lensing analysis taking into account 
the 3D geometry, that can be regarded as the true mass 
profile. Note that 3D mass profile for lensing has been re- 
covered via Abel inversion (assuming spherical symmetry) of 
the ID projected mass profile, while for X-ray-only data by 
assuming a spherical NFW. With respect to a joint analysis 
based on our triaxial modeling, a lensing-only analysis based 
on the standard spherical modeling predicts systematically 
higher masses by a factor 1.3—1.1 in the radial range out to 
400 kpc, moving from the center toward the SL boundary. 

An X-ray-only analysis based on the standard spherical 
modeling and strict HE predicts systematically lower masses 
by a factor of 1.5 — 1.2 in the radial range out to 1050 kpc 
compared to the true mass profile, moving from the cen- 
ter toward the X-ray boundary. Note that correcting the 



X-ray-only mass by £, e.g. by assuming a typical value from 
hydrodynamical numerical simulations, would lower this dis- 
crepancy by a factor ~ £. This confirms our insights about 
the role of the effects of geometry and non-thermal pres- 
sure on the physical properties. We also point out that the 
tree-dimensional mass discrepancy tends to get smal l er wit h 
increasing radii, as already pointed out by iGavazz 3 il2005h . 
This trend of the X-ray v.s. true mass discrepancy is a con- 
sequence of the fact that the IC gas is more spherical in the 
outer volumes than in the center, beside being altogether 
more spherical than the underlying DM halo: this is intu- 
itively understood because the potential represents the over- 
all average of the local density profile, and also because the 
gas pressure is isotropic unlike the possible anisotropic ve- 
locity ellipsoids for collisionless dark matter. 

We also point out that in IGavazz j (|2005T ) both the SL 
and X-ray 3D masses are systematically larger than the true 
3D mass distribution for clusters with prominent SL fea- 
tures, unlike in Figure [6] Indeed, he considered a toy model 
where the ICM and DM ellipsoids are azimuthally symmet- 
ric with the major axis oriented along the line of sight, the 
gas is wholly thermalized and it can be described as an 
isothermal /3-model, and all the physical parameters but the 
elongation along the line of sight are kept frozen. On the 
contrary, in the present analysis, owing to the full triaxial- 
ity, i.e. the two axis ratios and the three Euler angles, to 
non-thermal pressure support and to a physically-accurate 
description of the IC gas, we fitted this triaxial model to the 
data, allowing all the physical parameters to simultaneously 
change in the joint fit. 

We also compare the re s ults in our triaxial framework 
with those of iNewman et all l|201ll) , who adopted a triaxial 
DM distribution on the total mass distribution from lensing 
(assuming major axis of the DM halo along the line of sight) 
and with the radial profile of the stellar velocity disper- 
sion of the cD galaxy. They also accounted for hydrostatic 
mass estimates from X- ray data under th e assumption of 
spherical geometry from I Allen et al.l (|2008I ) . who used older 
Chandra calibrations with respect to those adopted in the 
present work, which lead to a higher spectral temperature 
(and therefore mass) profile by ~ 10%. They corrected the 
X-ray masses for violation of strict HE assumption by as- 
suming a typical value of £ ~ 10% from hydrodynamical 
numerical simulations. They assumed that the X-ray-only 
masses are recovering the true spherically-averaged masses, 
regardless of the halo shape. Whil e the results fro m hydro- 
dynamical numerical simulations (|Lau et al.ll2009T ) indicate 
that the bias involved in this assumption is just a few per- 
cent, this is strictly speaking correct for an average cluster, 
so this assumption likely merits some caution when we con- 
sider highly elongated clusters like Abell 383: as previously 
discussed, we found that the X-ray v.s. true 3D mass dis- 
crepancy is appreciable (see Figure [6}. 

Indeed, they found physical parameters dissimilar from 
ours: l/r)DM,a = l-97±g;?l, 7 = 0.59±g;lg and R s = 112±|J 
kpc. Note that we used a different definition of _R S , which 
refers to the majo r axes of the triaxial D M halo, so to make 
a comparison with lNewman et al.l (|201l[ ) we should multiply 
their _R S by ~ l/r}DU,a- still their _R S would be smaller than 
ours. 

The disagreement between the two works about physi- 
cal parameters, in particular with respect to the inner slope 
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Figure 6. 3D masses enclosed within a spherical aperture of ra- 
dius R from lensing-only data (long-dashed line with the 1-cr error 
dark gray shaded region) and from an X-ray-only analysis under 
the assumption of spherical geometry (short-dashed line with the 
1-cr error gray shaded region), and from a joint X-ray+lensing 
analysis taking into account the 3D geometry (solid line with the 
1-cr error gray shaded region). 



of the DM, suggests the follo wing considerat i ons: i ) there is 
the possibility that the fit of iNewman et all l)201ll ) is dom- 
inated by the high quality data for the circularized stellar 
velocity dispersion, therefore they might probe only the very 
central slope, whereas in the current paper we do probe the 
mass profile from 25 kpc out to ~1000 kpc; ii) for this clus- 
ter the gNFW may not be the correct densit y profile, which 
could have a slope that varies with radius (|Navarro et al.l 
|2004 ). though further data (e.g. stellar velocity dispersion) 
in a triaxial framework are needed to prove this: with this 
respect, both works might simply measure different part of 
the density profile; iii) the bias in the stellar velocity disper- 
sion resulting from proj ection effec ts (Dubinski 1998) might 
affect mass estimates of INewman et al" I (1201 ll ). 

A more thorough comparison between the two works is 
beyond the purpose of this paper, and it would involve a 
complex understanding of how their simplified assumptions 
e.g. on the X-ray data constraints, would affect the param- 
eter estimates. 



strated that accounting for the three-dimensional geometry 
and the non-thermal component of the gas allows us to mea- 
sure a central slope of the DM and concentration parameter 
in agreement with the theoretical expectation of the CDM 
scenario, dispelling the potential inconsistencies arisen in the 
literature between the predictions of the CDM scenario and 
the observations, providing further evidences that support 
the CDM paradigm. 

We also measured the amount of the non-thermal com- 
ponent of the gas (~ 10% of the total energy budget of the 
IC gas): this has important consequences for estimating the 
amount of energy injected into clusters from mergers, accre- 
tion of material or feedback from AGN. 

The increasing precision of observations allows test- 
ing the assumptions of spherical symmetry and HE. Since 
a relevant number of cosmological tests are today based 
on the knowledge of the mass and shape of galaxy clus- 
ters, it is important to better characterize their physical 
properties with realistic triaxial symmetries and account- 
ing for non-thermal pressure support. Galaxy clusters play 
an important role in the determination of cosmological pa - 
rameters such as the matter density (|Allen et al.l 1 20081 ), 
the amplit ude and slope of the density fluctuations power 
spectrum (Vocvodkin fc Vikhlininl 12004 ) . the Hubble con- 
stant illnagaki et al.lll995h. to probe the nature of the dark 
energy ( Albrecht fc Bernstein! 120071 ) and discriminate be- 



tween different cosmolog ical scenarios of structure formation 
l|Gastaldello et al.ll2007l ). It is therefore extremely important 
to provide a more general modeling in order to properly de- 
termine the three-dimensional cluster shape and mass, in 
order to pin down dark matter properties and to exploit 
clusters as probes of the cosmological parameters with a 
precision comparable to other methods (e.g. CMB, SN). 

The application of our method to a larger sample of 
clusters will allow to infer the desired physical parameters 
of galaxy clusters in a bias-free way, with important impli- 
cations on the use of galaxy clusters as precise cosmological 
tools. 
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6 SUMMARY AND CONCLUSIONS 

In this paper we employed a physical cluster model for 
Abell 383 with a triaxial mass distribution including sup- 
port from non-thermal pressure, proving that it is consistent 
with all the X-ray and SL observations and the predictions 
of CDM models. 

We managed to measure for the first time the full triaxi- 
ality of a galaxy cluster, i.e. the two axis ratios and the prin- 
cipal axis orientation, for both DM and ICM. We demon- 
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